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DATA FILTERING FOR THE NEW MINITRACK SYSTEM 

V. R. Simas 
D. E. Santarpia 

ABSTRACT 


Narrowband predetection filtering is required to avoid de- 
terioration in the SNR by the signal detector. The minimum pass- 
band, set by Doppler and other frequency uncertainties, is 7 KHz. 
Electronic post-detection filtering and data smoothing by polyno- 
mial fitting are equivalent, providing sufficient data samples are 
supplied. Since the data smoothing process exhibits an effect on 
the data equivalent to a filter having a passband of a fraction of a 
Hz, the electronic filter does not control the phase uncertainties. The 
proper data sampling rate is directly proportional to the j^st- 
detection bandwidth, B^. Because of this factor, should be 
made small avoiding excessive costs in relaying the data to GSFC 
and processing the data in the orbit determination sequence. Pas- 
sive filters for B^ cannot function without Introducing phase de- 
lays in the data very much greater than the desired accuracy of 
the system. For this reason, the tracking filter is selected which 
can follow, with a negligible error, the first, second and third de- 
rivatives of the signal phase resulting from a worst case space- 
craft pass of 100 miles altitude directly over the station. The loop 
bandwidth, B^, of this filter is a compromise between sampling 
rate and loop acquisition time. The best value is considered to be 
5 Hz which provides an acquisition time of 0.8 seconds and a 
sampling rate of 5 s/sec. 
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DATA FILTERING FOR THE NEW MINITRACK SYSTEM 


Introduction 

There are some interesting features in the signal filtering that occur in the 
Minitrack system. The three types of filtering that directly effect system per- 
formance are the pre-detection passband, the post-detection passband and 
finally the smoothing and polynomial fitting that is accomplished on the data 
by the computers at GSFC. 

The pre-detection filtering, although related to the other two, can be 
discussed independently. Its function is to reduce the noise on the signal that is 
subsequently demodulated by either an envelope or a square law detector. As 
the pre -detected signal -to-noise ratio (SNR) descends to unity and below for a 
non-coherent detector of this type, the SNR output of the detector becomes in- 
creasingly degraded. Thus, it is desired that the predetection filter be as 
narrow as possible. The width of this filter, however, must be sufficient to * 
encompass the Minitrack beacon signal in the presence of doppler and other 
frequency uncertainties. For the Minltrack parameters at 136 MHs, it can be 
shown that at least a 7 KHz passband is required. The present Mini track uti- 
lizes a 10 KHz linear phase filter. 

Post-detection electronic filtering and data smoothing are highly inter- 
related. From a theoretical viewpoint, it can be said that there is no necessity 
for electronic filtering because the smoothing operation which follows is in 
essence a filter with a very narrow passband, thereby precluding any effect 
the electronic filter has on the data. Practical considerations, however, 
warrant judicious selection of the electronic filtering characteristics for several 
reasons: (1) Phase readings are ambiguous in that they repeat every 360 degrees, 
thus, large phase fluctuations would introduce errors in the data in excess of 
the errors produced in quantizing the data, and would therefore require sub- 
stantially more complexity and expense in the data smoothing process, (2) 

Lightly filtered, thus noisy data, requires substantially higher sampling 
rates for signal processing, thereby increasing the expense of transmitting 
this data from the field stations to GSFC, (3) The cost of processing this lengthy 
stream of data by the orbit determination computers would be unjustifiably 
excessive. 

The extent of the electronic filtering is a compromise between the factors 
mentioned above and the errors and operational complications introduced by 
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the filter. The existing Minitrack utilizes a passive filter with a passband of 
10 Hz for the fine baseline data. For perigee passes, where the spacecraft 
altitude is approximately 100 miles, the rate of change in phase output of the 
interferometer, i, is on the order of 2.5 Hz or 900 degrees per second. As- 
suming the single tuned, 10 Hz wide, Minitrack Filter is tuned precisely to the 
difference frequency, a change in frequency, Af, of 2.5 Hz causes close to 25 
degrees of phase shift through this filter. This is a factor of about 75 times 
greater than the resolution of the system and produces an intolerable error 
even though the data is compensated for in the data processing sequences. If 
a passive filter were the only solution to this problem, then it would probably 
be better to widen it substantially beyond the existing 10 Hz passband even 
though the sampling rate would increase materially. A phase-lock-loop track- 
ing filter, because of its capability of providing narrowband filtering with 
vastly reduced phase shifts or time delay errors, appears to offer an ideal 
solution to this problem, assuming of course, that its characteristics are 
properly matched to the Minitrack signals. 


Phase- Lock- Loop Tracking Filter 

There is an apparent contradiction involved with the operation of tracking 
filters. Certain types of digital filters store data and can therefore look into 
the future as well as the past for their averaging functions thereby eliminating, 
for linear functions at least, time delay errors. These devices can, within 
limits, smooth out the noise fluctuations with no delay between the smoothed 
data and the applied data. Tracking filters achieve the same results in real 
time. From the vie^A^Doint of standard filter theory, this appears most incon- 
gruous, since conventional real time filtering is always involved with averaging 
past data so that the smoothed data is delayed by an amount corresponding to 
the filter weighting function. This is true even for data which changes linearly. 

An interpretation as to how the tracking filter accomplishes this feat can best 
be understood from an examination of Figure 1. The block diagram shows an 
idealized third order tracking filter in which the voltage controlled oscillator 
(VCO) is replaced by an integrator, a voltage controlled phase shifter (VC(^), 
and an oscillator. This manner of drawing the diagram portrays the integration 
that occurs due to the phase output of a VCO being the integral of the controlled 
frequency. We can now examine the operation of the circuit by applying waveforms 
to the loop 0 detector and deducing the relationship between the voltage analog 
of the phase variations of the applied signal and the VCO signal which is also 
the signal output from this filter. 

Consider applications of functions of increasing order to the loop. The 
first input is simply a phase offset. If the final integrator output, Vj(<^), is the 
proper voltage offset, the phase shift through the VC^will be sufficient to make 
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the tracking filter output (VC0 output) signal coincident and identical to the 
tracking filter signal input. To achieve this voltage offset, the input to the last 
integrator Vc (0)» will require zero voltage, and likewise for voltages and 
If they are not at the proper level, V^j (0) will be a voltage of higher 
order than that required to make the tracking filter output a replica of the 
signal input. The second row of waveforms shows a ramp change in phase 
applied to the tracking filter. In the same manner as before the waveform, 

Vj (0) must be a ramp change in voltage. This can be maintained by a voltage 
offset (<?!)). The voltages { 4 ) and (0) are zero to maintain a zero dif- 
ference between input and output. The third row of waveforms shows a quad- 
ratic change in applied phase. In the same manner as above, the output of the 
phase detector, (0), can now be zero. Finally for cubic changes in applied 
phase, the phase detector must produce a constant offset voltage, V^(0) whose 
magnitude is dependent upon the loop gain. Only for this last case is there a 
phase difference, hence a time offset, between the signal input and tracking 
filter output, this difference corresponding to the voltage offset at the output 
of the phase detector. 

It is important to realize that the deductions discussed above are 
predicated on steady state conditions only; that is, if a quadratic change in 
phase is assumed, then it has been changing a sufficient period, according to a 
quadratic relationship, for transient voltages in the loop to have asymptotically 
decayed to zero or to a steady state value. It is also important to note that a 
signal input can be highly filtered in this tracking filter with no phase or time 
delay, providing changes in applied phase (or frequency) occur sufficiently 
smoothly that third order terms and higher are negligible. 

As further evidence, let us view the effect that an artificial loop time delay 
would have on the difference between output and input. Consider, for example, 
the second row of waveforms where the applied signal contains a ramp change 
in phase. If a loop time delay were introduced after the last integrator it would 
appear at first glance that a change in phase would occur between the signal 
input and tracking loop output since the ramp voltage fed back to the VCO would 
be delayed. The fact is, however, the voltage, (0) , after the transient has 
decayed, will vary (after the delay) to reduce the phase difference between the 
signal input and the tracking filter output to zero. Thus, the circuit even in the 
case of artificially induced loop delays produces control voltages which can be 
considered to anticipate the proper value occuring in the future. This set of cir 
cumstances becomes plausible when considering the nature of the loop input 
signals. Since we are considering only well-defined waveforms in steady 
state, the loop has a priori knowledge and will track, without phase lags, wave- 
forms which can be described by equations whose order is one less than the 
number of loop integrators. For phase fluctuations of no higher order than 
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quidratic, our idealized filter, will anticipate perfectly the future eequence of 
events and will compensate accordingly* 

Because of the unique properties of tracking filters, as discussed above, 
such a device appears to be most appropriate for the Minitrack post-detection 
filter where it is desired that namw band filtering be accomplished with 
negligible phase delay* 

Another source of phase error which a trucking filter must accommodate 
is the dynamic phase error introduced by the time derivatives of the phase 
output of the Minitrack interferometer system, These time derivatives are 
functions of spacecraft-earth station geometry and spacecraft dynamics* It is 
important, therefore, to determine the maximum total phase error possible 
due to a severe satellite tracking situation and the tracking filter bandwidth 
required to maintain a phase error tolerance of, say, one count or ,36^. 


Tracking Dynamic Phase Errors 

Consider a worst case dynamic tracking situation in which a spacecraft 
orbiting at 5 mi/sec passes overhead in a 100 mile high circular orbit. Refer- 
ring to Figure 2: 

r = earth radius = 4,000 miles 
h - orbit height = 100 miles 

R = r + h 

= fipacecraft tangential velocity (5 miles/sec) = Ro) = -Ra 
a, d, p - are geometric angles determined by the spacecraft, 
station and earth in radians. 

From the geometry of Figure 1: 


and 


but 


y8 = 90 a - ^ 
yS = - a - ^ 



( 1 ) 


( 2 ) 
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Substituting equation (2) into (1): 


V, . 
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also 




From the interferometer relationship:^^ ^ 


cos ^ ^ 
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where; 


B == Minitrack interferometer baseline of 50 wavelengths 

- Phase output of the minitrack interferometer system in cycles. 

but from Figure 1; 


cos 0 


R . . 

— Sin p 
r 


therefore, 


4 > = 


BR 


Sin 


ft 


differentiating once and substituting equation 3 for ft: 


i — cos — fi cos /3 

T r 


( 7 ) 


taking the second derivative: 


$ = s in /3 (- V. + R^) - . — cos /3 (8) 

r ' r 


Substitution of equations (4), (5) and (6) into (8) gives the relationship of 
the phase acceleration, to spacecraft-earth tangential angle, 6, A solution to 
equation (8) was obtained by computer for 0 ^ 6 ^ n in .01 radian increments. 

In addition, computer differencing was used to determine the rate of change 
of phase acceleration, . A plot of $ and as a function of spacecraft- 
earth tangential angle, 0, is shown in Figure 3. 
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Determination of Phaae-i:x)ck-Loop Tracking Bandwidth 


In a system such as the 136 MHz Minitrack, the steady state phase error, 
in radians, introduced in the tracking loop due to « and 4>' is described by 
Reedy and Mallinkrodt^ to be: 


e(2"** order) 



2/2 V 




oz 


( 6 ) 


e(3"* order) = 


%7T •!* 


n 


(7) 


where 'i, 4>' are time derivatives of the interferometer phase (in 

cycles/sec^ and cycles/sec^ respectively) at the input to the loop phase detector, 
and is the "natural frequency" of the loop and is equated to the one sided 
noise bandwidth of the loop by^*^ : 


order) 



(for damping factor .707) 


order). = 


.73 


(Mallinckrodt Loop) 


substituting this relationship for a , we obtain: 


e(2"‘‘ order) = 


B* Bl 


( 8 ) 


eCd'** order) = 


9.78 -i- 


(9) 


Maximum total steady state phase error for a second order tracking loop 
occurs for a spacecraft at ±27® off zenith. At this spatial position « is maxi- 
mum and ^ = 0. For a third order loop, the maximum phase error occurs at 
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zenith where <t> is maximum; hence equations (8) and (9) reduce to: 


6 (2"*^ order) = 

max ' ' 


1.77 




( 10 ) 


^max order) 


9.78 

<P 

■» " 
Bl 


( 11 ) 


solving for we obtain: 


0^(2"** order) 


1.77 

€ 

max 



(12) 


^(3'** order) 



(13) 


• • • • • 

The maximum values of <l> and 4> from Figure 3 are .0987 cycles/sec^ and 
.0167 cycles/sec^ respectively. Inserting these values in the above equations 
with set to 0.36 degrees (1 count) results in: 

(2nd order) = 5 Hz. 

Bj^ (3rd order) = 3 Hz. 

Having set a lower bound for selecting acceptable phase lock loop tracking 
filter bandwidths, it is necessary to examine the effect B^ has on data acc|ui- 
sition time and the subsequent smoothing operation performed by the GSFC 
orbit determination computer. 


Acquisition Time and Data Smoothing 

Acquisition time is a most important parameter in the selection of accept- 
able phase-lock-tracking filters. Since a 100 mile orbiting spacecraft is in the 
main lobe of the minitrack system, only 4 or 5 seconds the acquisition time 
should be on the order of 1 second. Assuming the interferometer system 
begins operation as a spacecraft just enters the main lobe, approximately 5® off 
zenith, the following phase error-producing interferometer phase derivatives 
exist: 
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= - .225 rad , 

♦ = 15.07 

sec 

♦ = .1858 

sec^ 

*4)' = -.0992 -L®i 
sec^ 


As the satellite approaches zenith^ 4> and *4»* are increasing to a maximum 
value and $ is decreasing to zero; however, at commencement of system opera- 
tion, transient phase errors in the tracking filter are produced which require a 
finite time to decay to the desired steady state resolution of 0.36^. On the 
basis of the transient res ponse of a 5 Hz second order phase lock loop whose 
dampening factor is /2/2, approximately 0.8 second is required to reduce 
tracking loop transient phase errors to .36^. 

The transient response of a third order loop, given the same initial phase 
error producing parameters given above, is greater than four seconds^ which 
is a large portion of the main lobe tracking time. 

It is apparent that a 5 Hz second order loop will adequately meet the 
criteria set for acceptable acquisition time, while a 3 Hz third order loop will 
not. The third order bandwidth can be widened to reduce the loop transient 
error response time to less than a second at the expense of S/N ratio improve- 
ment and increased sampling rate. About 20 Hz would be required. It is also 
possible to attempt acquisition earlier in the pass, that is, prior to the main 
lobe beam width. If acquisition is attempted prior to the main lobe null, 
questions arise as to the magnitude of the phase errors introduced due to per- 
turbations experienced in passing through the lobe nulls. Would these phase 
errors be of sufficient magnitude to produce a loss of lock? Would they be of 
sufficient magnitude to produce longer transient phase errors ? For these 
reasons the 2nd order loop with a 5 Hz bandwidth is the most appropriate for 
our application. 

The smoothing program performed at GSFC takes the total data message 
consisting of N data points and by a polynomial fitting process determines 
the spacecraft trajectory arc over which the data was observed. The uncertainty 
of the smoothed data sent to GSFC can be reduced by a factor very nearly equal 
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to (where N is the number of samples) providing the sampling rate 

is sufficiently low that the noise on these samples is uncorrelated. The noise 
starts to become correlated when the sampling rate is on the same order as the 
post detection bandwidth, \ . Increasing the sampling rate beyond this point 
results in a negligible improvement of the processed data. 


Trade-Offs and Conclusions ; 

It is apparent from the preceeding that the selection of a tracking filter 
which will meet the demands of the Minitiuck System is made on the basis of a 
compromise between sampling rate, f^ , and acquisition time. For economic 
reasons, sampling rate should be as low as possible while still maintaining the 
sampling criteria for processing data. As discussed earlier, increasing the 
sampling rate beyond does not improve the processed data; therefore, we 
recognize the necessity for selecting a tracking filter on the basis of band- 
width since the narrower the bandwidth required to maintain a desired 
Minitrack resolution, the lower the required sampling rate for processing the 
data; however, it must be remembered that the smaller the bandwidth, the 
longer transient errors exist. 

As previously showsi, a second order tracking filter of 5 Hz or greater 
and a third order tracking filter of 3 Hz or greater will adequately track 
steady-state phase errors while maintaining a residual error of .36®. Because 
of uncertainties involved in attempting acquisition prior to the main lobe null, 
it is advisable to count on acquiring only after this null point is passed. A 3 Hz 
third order tracking filter will not perform with this constraint because of the 
duration of its transient response. Increasing the bandwidth to 20 Hz^®^ alleviates 
the problem but at the expense of added noise and increased sampling rate. A 
5 Hz second order tracking filter will provide the necessary tracking accuracy 
and acquisition time without going to the wider third-order loop and higher 
sampling frequency. 

In view of the above considerations, appropriate values for a tracking 
filter bandwidth and sampling rate are considered to be 5 Hz and 5 samples per 
second providing the type of tracking loop utilized is second order. 
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